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Abstract 



A fractal oscillatority of solutions of second-order differential equations near 
infinity is measured by oscillatory and phase dimensions. The phase di- 
mension is defined as a box dimension of the trajectory (x, x) in M 2 of a 
solution x = x(t), assuming that (x,x) is a spiral converging to the ori- 
gin. In this work, we study the phase dimension of the class of second-order 
nonautonomous differential equations with oscillatory solutions including the 
Bessel equation. We prove that the phase dimension of Bessel functions is 
equal to 4/3, and that the corresponding trajectory is a wavy spiral, exhibit- 
ing an interesting behavior. The phase dimension of a generalization of the 
CN ■ Bessel equation has been also computed. 

l> 
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1. Introduction and motivation 



The fractal oscillatority of solutions of different types od second-order 
linear differential equations has been recently considered by Kwong, Pasic, 
Tanaka and Wong. The Euler type equation has been studied in |7( and J8j, 
the Hart man- Wintner type equations in [5], half linear equations in 10 1, and 



finally the Bessel equation in |9j . In all of this work the fractal oscillatority 
is considered in the sense of the oscillatory dimension, at first introduced in 
Pasic, Zubrinic and Zupanovic (ll|. The oscillatory dimension of a solution 
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x(t) is defined as the box dimension of a graph of the function X(t) : = x(-) 
near r = 0. 

On the other hand, the fractal properties of spiral trajectories of dynami- 
cal systems in the phase plane have been studied by Zubrinic and Zupanovic, 



see e.g. [13[ and [15l . From their work the co ncep t of the phase dimension 



has arisen and has finally been introduced in [11]. They adapted standard 
idea of phase plane analysis to fractal analysis of solutions of second-order 
nonlinear autonomous differential equations. 

These ideas motivated us to study a fractal connection between the os- 
cillatory and phase dimensions for a class of oscillatory functions, see [3]. In 
that study, we discovered a specific type of spirals with a nondecreasing ra- 
dius function, related to chirp-like solutions of a class of equations considered 
by Kwong, Pasic, Tanaka and Wong, which we call the wavy spirals. 

A model for chirp-like behavior of solutions developed in [3[ could not 
handle the specific behavior of Bessel functions, so to determine the phase 
dimension of Bessel functions we generalize the technique used in [3[ . Notice 
that the oscillatory dimension of Bessel functions was already considered in 

0- 

In this work we actually consider some generalization of the Bessel equa- 
tion that is motivated by the generalization introduced in [9j, for which we 
determine the phase dimension in our main result, Theorem [TJ The standard 
Bessel equation now becomes a specific case of this generalization. In order 
to p rove Theorem [T], we first obtain a new version of some theorems from 



13 



We discover an interesting property related to wavy spirals. By varying 
values of parameters in our generalized Bessel equation, we go from spirals 
with no waves to spirals with "big" waves. 

2. Definitions 

We first introduce some definitions and notation. For A C M. N bounded 
we define the e -neighborhood of A by A e := {y e M. N : d(y, A) < e}. By the 
lower s- dimensional Minkowski content of A, s > we mean 



v 



Ml{A) := liminf 



\A 



: I 



e^O e N ~ 



and analogously for the upper s- dimensional Minkowski content A4* S (A). 
Now we can introduce the lower and upper box dimensions of A by 

dinipA := inf{s > : Ml(A) = 0} 



and analogously dim^A := inf{s > : M.* S (A) = 0}. If these two values 
coincide, we call it simply the box dimension of A, and denot e by dini p A . 



For more details on these definitions see e.g. Falconer [2j, [13[ and [14 . 

Assume now that x is of class C 1 and t > 0. We say that x is a phase 
oscillatory function if the following condition holds: the set T = {(x(t), x(t)) : 
t e [to, oo)} in the plane is a spiral converging to the origin. 

By the spiral here we mean the graph of a function r = f(<p), <^ > <^i > 0, 
in polar coordinates, where 

/ : [if i, oo) — > (0, oo) is such that f(ip) — > as ip — > oo, 

/ is radially decreasing (i.e., for any fixed ip > <fi (1) 

the function N 3 k >->■ /(</? + 2&7r) is decreasing). 

Depending on the context, by the spiral here we also mean the graph of a 
function r = #(<£>), (p < <fi < 0, in polar coordinates, where for h((p) = g(—(p), 
V</5 > \Pi\i the graph of a function r = /i(y?), v 9 ^ Ivil > 0? given in polar 
coordinates, satisfies ([T]). It is easy to see that the spiral given by function g 
is a mirror image of the spiral given by function h, regarding x-axis. We also 
say that a graph of function r = f(<p), tp > (fi > 0, in polar coordinates, is 
a spiral near the origin if there exists (p% > <p\ such that a graph of function 
r = f(<p)i V 9 ^ ^2 is a spiral. 

The phase dimension dim p h(x) of function x : [to, 00) — ¥ M. of class 
C 1 is defined as the box dimension of the corresponding planar curve T = 
{(x(t),x(t)):te [t , 00)}. 

We use a result for the box dimension of spiral T defined by r = <p~ a , 
V 9 ^ V^o > 0, dims T = 2/(1 + a) when < a < 1, see Tricot [12,, p. 121] and 
some generalizations from [131 ] . 

The phase dimension is a fractal dimension^ introduced in the study of 



chirp-like solutions of second order ODEs, see [ill ]. Fractal dimensions are a 



well known tool in study of dynamics, see [16 . 

For two real functions f(t) and g(t) of a real variable we write f(t) ~ g(t) 
as t — ?• (as t —7- 00) if there exist positive constants C and D such that 
C f(t) < g(t) < D f(t), for all t sufficiently large. For example, for a function 



F : U ->■ V with U,V CR 2 , V = F(U), the condition \F( Xl ) - F(x 2 )\ ^ 
|xi — X2I means that / is a bi-Lipschitz mapping, i.e., both F and F~ x are 
Lipschitzian. 

We write f{t) ~ g(£) if f(t)/g(t) — >• 1 as t — >• 00. Also, if A; is a fixed 
positive integer, for two functions / and g of class C k we write, 

f(t) ~ fe gf(t) as t ->• 00, 

if /&)(*) ~ 00) (f) as t -> 00 for all j = 0, 1, ..., k. 



i — a 



For example, - — X — ~3 t a as t — > 00, for a G (0, 1). 

We write f(t) — 0(g(t)) as t — > 00 if there exists a positive constant C 
such that |/(t)| < C|gr(i)| for all t sufficiently large. Similary, we write f(t) = 
o(g(t)) as t — > 00 if for every positive constant e it holds \f(t)\ < s\g(t)\ for 
all t sufficiently large. 

3. Phase dimension of Bessel functions 

From the point of view of a fractal geometry we study a spiral generated 
by Bessel functions and generalized Bessel functions. 

The Bessel equation of order v, widely known in literature, see e.g. 
p. 98], is the linear second-order ordinary differential equation given by 

t 2 x"(t) + tx\t) + (t 2 + u 2 )x(t) = 0, (2) 

with one parameter v G M.. It has two linearly independent solutions, known 
as Bessel functions of the first and second kind of order u, designated J u and 
Y u , respectively. 

Pasic and Tanaka in [9J] studied the oscillatory dimension of the Bessel 
equation and one specific generalization of this equation that they introduced. 
Motivated by their generalization we introduce this class of equations 

t 2 x"(t) + t(2 - /i)x'(t) + (t 2 + v 2 )x(t) = 0, (3) 

where i/,/i6l are parameters. By setting fi = 1 we get the standard Bessel 
equation (T5]). Two linearly independent solutions to equation ([3]) define func- 
tions J v ^ and Y U41 , which we call generalized Bessel functions. By solving (EJ), 
we see that J u ^ and Y v>li can be written in terms of J„ and Yp, respectively, 



l\ 2 



where v = v/(^) + v 2 ■ We get 



M -l 



JvAt) = t—Jv(t) 



M-1. 



^(i) = t—Y v (t). 

As /„(£) = O U~*\ and y„(i) = O (*"*), for every z/ £ R, it follows that 

Jv,fj,(t) = O [t^~ ) and Y V)li (t) = O [t^ ), for every u, fj, £ R. Now, assume 
that To > 0. It is easy to see the following: 

• If \x > 2 then J u ^{t) and y^ iM (t) are unbounded for £ £ [r , oo). Trivially, 
Jv,fi(t) and l^ >A j(t) are not phase oscillatory functions. 

• If /j, — 2 then J v ^{t) and l^ )M (t) are bounded for £ £ [r ,oo), but 
lim^oo J u ^{t) and lim^oo Y VjfM (t) do not exist, so J u ^(t) and ^(t) are 
not phase oscillatory functions. 

• If jjl < then J u ,fi{t) and 5^, M (£) are phase oscillatory functions, but 
with a trivial phase dimension, dim p h [J V A = dim p h [J^n) = 1, for 

every j/6l. 

The remaining case fi £ (0, 2), which will produce phase oscillatory functions 
with a nontrivial phase dimension, we describe in Theorem [TJ 

Theorem 1. (Generalized Bessel spiral) Let fi £ (0,2), v £ R and r > 0. 

Let x(t) = J Vt n(t), for t £ [r , oo), and let continuous function ip{t) be given 
bytmtp(t) = |||, fort £ [r ,oo). 

Then the planar curve T = {(x(t),x(t)) : t £ [r ,oo)} is a spiral near 
the origin r = g(<f), f £ (— oo, — </>o], satisfying g(ip) ~ |y9|~^~, in polar 
coordinates. Phase dimension 

4 
dim p / l (x) = dune T 



A-n 



Spiral V is a wavy spiral if and only if v ^ 0. Analogous claim is valid for 
x(t) = Y„ tli (t), t £ [t ,oo). 

See Definition [2] below for the wavy spiral. From Theorem [1] directly 
follows our result about Bessel functions. 





(a) Curve Ti, dime r\ = ^ 



(b) Graph of function x\(t) 



Figure 1: Curve Ti = {(xi(t) , Xi(t)) : t G [10, oo)} and a graph of function 
x\{t) = J 5 ,o.2(t). 





(a) Curve 1^, dims r 2 = | 



(b) Graph of function X2 (t) 



Figure 2: Curve T 2 = {(x2{t) , X2{t)) : t G [10, 00)} and a graph of function 
x 2 {t) = J 5 ,i(t). 





(a) Curve T 3 , dims r 3 = fy 



(b) Graph of function X3(t) 



Figure 3: Curve T 3 = {(x 3 (t) , x 3 (t)) : £ G [10, 00)} and a graph of function 
x 3 (t) = J 5 ,i.sW- 
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Corollary 1. (Phase dimension of Bessel functions) Phase dimension of 
Bessel functions is equal to 4/3. 

Proof. Follows directly from Theorem [T] by taking /j, = 1. □ 



Remark 1. It is interesting to notice that value 4/3 also appears as the box 
dimension of a trajectory near the origin of system 



r(r 2 + a ), a E 



<P 



in polar coordinates, by taking a® = 0, see [13j. This system is the normal 
form of the Hopf bifurcation, which is a well known bifurcation of 1-parameter 
families of vector fields. 

Another place where value 4/3 appears is as the box dimension of the 
famous Euler spiral, or the Cornu spiral, also called the clothoid, defined 
parametricaly by 

x(t) = / cos(s 2 ) ds, 
Jo 

y(t) = / sin(s 2 ) ds, where t E M, 
Jo 

near t = ±oo, see J4j. The Euler spiral appears in various branches of 
mathematics, physics and engineering. Some applications are in optimal 
control theory, in the description of diffraction phenomena in optics and in 
robotics for pathfinding algorithms. 

In Theorem [1] we have a spiral generated by generalized Bessel functions, 



see Figures dHSl To prove Theorem [TJ we need a new version of 13|, Theorem 



5] with weaker assumptions on the spiral than in 13 



Theorem 2. (Dimension of a piecewise smooth nonincreasing spiral) Let 
f : [(pi,oo) — > (0, oo) be a nonincreasing and radially decreasing function, 
also a continuous and piecewise continuously differentiable. We assume that 
the number of smooth pieces of f in [(fij^fx) is finite, for any Tp x > tp x - 
Assume that there exist positive constants m, m, a and M such that for all 

rrup~ a < f(ip) < fnip~ a , 



a^- a - 1 </(^)-/(v + 27r), 
and for all ip where f(ip) is differentiable, 

\f^)\<M^ a -\ 

Let T be the graph of r = f(ip) in polar coordinates. If a £ (0, 1) then 

2 



dims T 



1 + a 



We also need the following Lemma [U that is a generalization of 13 . 
Lemma 1] dealing with smooth spirals. 

Lemma 1. (Excision property for piecewise smooth curves) Let T be the 
image of a continuous and piecewise continuously differentiable function h : 
[if i, oo) — > M? (piecewise in the sense of Theorem^. Assume that dim B r > 
1, r\ := h((Tp l7 oo)), for some fixed Tp x > ip 1 , and h{\ip\ : Tp^\) f]Ti =0. Then 



dim^r! = dim s r, dim^T! = dim B r. 

Notice that T from Lemma[T]is a piecewise smooth curve in M 2 . Theorem[2] 
and Lemma [T] are introduced and proved in [3J]. 

Actually, we need a version of Theorem [2] with even weaker assumptions 
on associated spiral. So, using Theorem |2] and Lemma (TJ here we prove 
Theorem [3] that deals with a spiral T' described by r = f{ip), where / is 
increasing on some parts, see Definitions [1] and [2j We call this new property of 
r" spiral waviness, which naturally arises in spirals generated by generalized 
Bessel functions in Theorem [TJ 

Definition 1. Let r : [to, oo) — > (0, oo) be a C 1 function. Assume that 
7"'(^o) < 0. We say that r = r(t) is a wavy function if the sequence (t n ) 
defined inductively by: 

t 2k+1 := inf {t : t > t 2k , r'(t) > 0}, ke N , 

hk+2 ■= inf {t : t > t 2k +i, r(t) = r(t 2k +i)}, k £ N , 



is well-defined, and satisfies the waviness condition: 

' (i) The sequence (t n ) is increasing and t n — > oo as 

n — > oo. 



< 



(ii) There exists e > 0, such that for all k G No holds 

^2fc+l — hk > £ - \ ' 



m 



For all fc sufficiently large it holds 

osc r(t)=o(t- k a - 1 1 ),ae (0,1), 



where oscr(t) = maxr(t) — minr(t). 
tei t&i tei 



Notice that min r(t) = r(t2k+i)- Condition (i) means that the 

iG[*2fc+l)*2fc+2] 

property of waviness of r = r((p) is global on the whole domain. Condi- 
tion (ii) is connected to an assumption of Lemma [2j Condition (Hi) is a con- 
dition on a decay rate on the sequence of oscillations of r on Ik = [^+1,^+2], 
for k sufficiently large. Also, notice that condition r'(t ) < assures that t\ 
is well-defined. 

Remark 2. Conditions (i) and (ii) in the waviness condition (jj]) are not 
entirely independent. From (ii) and if (t n ) is increasing follows that t n — > 00 
as n — > 00, but from (i) does not follow (ii). So, condition (ii) plus (t n ) 
increasing is stronger than condition (i). 

Definition 2. Let a spiral T', given in polar coordinates by r = f(ip), where 
/ is a given function. If there exists increasing or decreasing function of class 
C 1 , ip = (f(t), such that r(t) = f(<p(t)) is a wavy function, then we say T' is 
a wavy spiral. 

Now, using Theorem [2] and Lemma [1] we prove the following Theorem [3j 

Theorem 3. (Box dimension of a wavy spiral) Let to > and assume r : 
[t , 00) — )■ (0, 00) is a wavy function. Assume that ip : [t , 00) — > [(p , 00) is 
an increasing function of class C 1 such that tp(to) = <fo > and there exists 
^o6ffi such that 

\(<f(t) - (po) — (t-t )\ -)• as t — ¥ 00. (5) 



Let f : [ipo, oo) — > (0, oo) be defined by f{(f{t)) = r(t). Assume that V is a 
spiral defined in polar coordinates by r = f{(f), satisfying (JTJ). Let a G (0, 1) 
is the same value as in (jl])(iii) for wavy function r, and assume e' is such 
that < e' < e, where e is defined by (jl])(ii) for wavy function r. Assume 
that there exist positive constants m, m, a' and M such that for all ip > <p 0! 

mp' a < f(<p) < rn^ a , (6) 

\n<p)\<M<p~-\ (?) 

and for all /Sip, such that 6 < Aip < 2tt + 6, there holds 

oV""- 1 </(¥>)-/(¥> +A¥>), (8) 

where 9 := min {e', it}. 

Then V is a wavy spiral and 

diniB T' 



1 + a 

Proof. 

We will construct a new spiral T that is nonincreasing and close to V 
in the sense that there exists Bilipschitz map F : V — > F (ie, there exist 
constants K\ e (0, 1) and K 2 > 1 such that 

K x d(X',Y') < d(F(X'),F(Y')) < K 2 d(X',Y'), 

for every X', Y' G T', where d(X', Y') is Euclidian distance between points 
X' and Y'). Also, we will construct this map F. We know from [2|, p. 44] 
that then dmie V = dim^T. Because of Lemma [1] we have to construct T 
and establish Bilipschitz map F only for cpo and equivalently to sufficiently 
large. This means we can increase value to to t 2 k and take <po = <p(t2k ) f° r 
any /co G N, where sequence (t„) is from Definition [TJ 

Let T be defined by r = f(<p), in polar coordinates. Informally, T is 
composed of alternating radially nonincreasing parts of spiral r" (with respect 
to parameter t) and circular arcs that substitute the parts where T' radially 
increases and then decreases until it intersects the associated arc. We see 
that T is radially nonincreasing. 

Formally, we define T by parts as follows: 

r|[t 2fc ,t 2fc+1 ] := r'| [t2fe><2fc+l] , for every A; G N , 
f \[t 2k+1 ,t 2k+2 ] ■= {{r{t),<p{t)) : r(t) := r(t 2k +i), t G [t 2k+1 , t 2k+2 ]} , 
for every k G No, 
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in polar coordinates. 

As r is a wavy function (see Definition [IJ and T' is a spiral, satisfying 
([[]), we see that T is well denned and is also a spiral, satisfying ([[]). Notice 
that f(tp(t)) = r{t) and r(t) < r(t), for all t G [t , oo). 

Assume X" G T' , which means T" is given by (r(£jv); v(£f')) f° r some value 
£t', in polar coordinates. Now we define 

F(T'):=(r(i r 0^(iT0), 

in polar coordinates. We see that map F radially maps points from spiral r" 
to spiral T. 

We will use f(T') and f(T) to denote radius, in polar coordinates, of 
points T' G r' and T G T, respectively. With <p T we denote argument, in 
polar coordinates, corresponding to points T" G V or T = F(T') G I\ With 
O we denote the origin of the phase plane. 

Now, we prove that map F is a Bilipschitz map. Notice this is equivalent 
to the claim that F' 1 : T — > r" is a Bilipschitz map. 

We have five cases depending on the relative position of points X,Y G T. 
For each case we will prove that F _1 is a Bilipschitz map. 

Case A. Assume point X G r|r t2fc)t 2k+1 i, for some k G N , and point Y G 
r|[t 2l ,t 2!+1 ], for some Z G N . Map F" 1 is identity map on f\^ 2kMk+l] U? | [42; , 42;+l] 
so F~ l is, trivially, a Bilipschitz map on T\[t 2k ,t 2k+1 ][jT\[t 2l ,t 2l+1 ], for every 
Me N . 

Case B. Assume points X, Y G r|r t ;t2 i, for some fc G No, and 
Iv 9 ^ — ¥>y\ < 0. Without loss of generality we can assume that ipx < <£y- Let 
us denote with X' = F _1 (X), Y' = F _1 (y), and with X x (X 2 ) the point on 
spiral T corresponding to parameter t 2 k+i (^2^+2) and angle ^2^+1 : = ffek+i) 
(f2k+2 '■= ^p(t2k+2)), see Figure HI As ip is increasing, we see that (p2k+i < 

i P2k+2- 

Using elementary geometry regarding isosceles triangle XYO and triangle 
X'Y'O (see Figure HJ) we conclude that 

d(X,Y) < d(X',Y') = d{F- 1 {X),F-\Y)). 



Now, let us assume that /(^y) > f{fx) and let point Z on the line OF' 
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Figure 4: Geometry in the proof of Theorem [3], Case B. 



be such that d(0, Z) = d(0,X r ). We see that 

d{F-\X),F- 1 {Y)) = d{X',Y')<d{X',Z) + d{Z,Y') 

< \ Vx -ip Y \-d(p : X') + \f(X')-f(Y') 

= \(px-<PY\-f(<Px) + \f(<Px)-f(<PY)\ 

= \ip x -ip Y \- f(<Px) + \f(0\\ ( Px -<Py\, 



< 



l^x - <£y| 



/(¥*) + sup |/'(0| 

?G[l/32fc+l,V 5 2fc+2] 



l^x - fvlfi^f 



2k+l) ' 



/(</?*) + sup \f (0| 

££[V2fc+l>¥>2fc+2] 



/(^2fc+l) 

Analogously, if f(<fy) < f{fx) we conclude that 

f(<py)+ SUp |f (0| 



d(F-\X),F-\Y)) < \<px- <pr\f (<P2k+i) 



/(V2fc+l) 
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As 9 < it, using Proposition [TJ we see that 

\<Px ~ <PY\f(<P2k+i) < ^d{X, Y). 

As function / satisfies conditions flS} and (J7J), and as we can take (p 
sufficiently large such that Mip^+i < 7 rn\ l f2k+i\~ a , we compute 



f(<fx) + sup If (0| mpx a + sup (M£ 

£,£[<P2k+l,<P2k+2] , i&[V2k+1^2k+2\ 



-a-l\ 



f(<P2k+i) ' mP2k+i 



< ^2W + ^2fc+l 1 < m f2k+l + m ^2k+l 
lW>2k+l " W+Pzk+l 













< 


m 
2— 














m 


Analogously, 


we 


get 


fifPv) 


+ 




sup 



l/'(0l 

?S[V2fe+li¥'2fe+2] ,, n fn 



< 2- 



/(^2*+i) 2a' 

Finally, we conclude that 

TV TH 777 

^(X), F-!(y)) < -d(X, y) • 2- = 7T-d(X, Y), 

2 mm 

so F _1 is a Bilipschitz map on r|[t 2fc+1 ,t 2fc+2 ], for every k G No. 

Case CI. Assume point X G r| [t2fc+2i<2fc+3 ] and point Y G f% 2ifc+1 ,t 2fc+2 ], 
for some k G No, and \<px — 9?y| < 0. Let us denote with Z G T the point 
corresponding to parameter t 2 k+2 and angle y?z- Notice that (px > <fiz > <Py- 

Let us denote with X' = F _1 (X), V = F" 1 ^) and Z' = F" 1 ^), see 
Figured Notice that X = X', Z = Z', which gives d(X,Z) = d(X',Z'). 
Now, 

f(X) = d(0,X) = d(0,X') = f(X f ), 

f(Z) = d(0,Z) = d(0,Z') = f(Z'), 

d(0,Y) = f(Y)<f(Y') = d(0,Y'), 

and notice that 

d(0,Y) = f(Y) = f(Z) = d(0,Z), 

13 




From 



Figure 5: Geometry in the proof of Theorem [3j Case CI. 

f(Z') = f(Z) = f(Y) < f{Y'). 

d(Q,X) = f(X) < f(Z) = f(Y) = d(0,Y) 



easily follows (see Figure |5]) that the angle between line XY and line YY' in 
triangle XYY' is greater than or equal to tt/2, so using elementary geometry 
we see that 

d(X,Y) < d{X\Y') = d(F- l (X),F-\Y)). 

On the other hand, from Case B, it follows d(Z',Y') < ir^d(Z,Y). 
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Now, we compute 

d{F-\X), F-\Y)) = d(X', Y') < d(X', Z') + d(Z', Y') 

777 777 

< d(X, Z) + Tr—d(Z, Y) < it- [d(X, Z) + d(Z, Y)] 

m m 

rn 

< it— [\<p x - Vz\ max{d(0, X), d(0, Z)} + (d(0, Z) - d(0, X)) 
+ \<p z -(pY\Tnax{d(p,Z),d(p,Y)}] 



m 



7T- 



m l 
m 



(\<Px - <Pz\ + \<pz ~ <Pv\)f(Z) + (f(Z>) - f(X')) 



7T- [\<p x - V y\f(Z') + (f(<p z ) - f(<p x ))] 



< 7T- 



m 
m 
m 

m , 



\<fx-<fY\f(<fz) + \<fx-(fz\ sup |f(OI 

££[<Pz,<Px] 



< k—\<Px-<Py\ 
m 



f(<py)+ SUp |f (0| 



n— ■ \<p x - <Py\J{Vx) ■ 77 ; 

m fvpx, 



/{<&)+ sup If (01 

ffar) 



Let point Y" E OY be such that d(0,Y") = d(0,X). Notice that 
d(0, Y") < d(0, Y). As before, using elementary geometry regarding isosce- 
les triangle XOY" and triangle XOY (see Figure[5]) we conclude that d(X, Y") < 
d(X,Y). Using Proposition [1] on isosceles triangle XOY", we get that 

\<Px-tPY\f{<Px) < \d{X,Y") < ^d(X,Y). 
We can take (p sufficiently large such that, using <px > ¥>y > y?o 5 we get 
f{<p Y ) , m<fy c 



/(? 



x, 



< 



< 



mip x m(ip Y + \<fix ~ ^y\T 



rrup Y 



m(<p Y + 0) a _ rnipyil + 0(^ 1 )) 



<2™. 
m 



rmpy rrupy 

Analogously as in Case B, we can take <po sufficiently large such that 

f(<py)+ SUp |f (01 

Z&[<pywx] < m 
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Figure 6: Geometry in the proof of Theorem [31 Case D. 



Finally, we see that 



d(F-\X), F-\Y)) < tt™ ■ U{X, Y)-2™. 2™ 

m 2 mm 



2it 2 



d(X,Y). 



Case C2. Assume point X £ F\[ t2k+ljt2k+2 ] and point Y £ F\[ t2k>t2k+1 ], for 
some k £ N , and \tpx ~ Vy I < 0- Proof is analogous to proof of Case CI. 

Conclusion of Cases A— C2. Using (@J(ii) and ©, as t 2 k+i—t2k > £ > 
for all k £ N , it follows that, for t and equivalently <p sufficiently large, 
cases A to C2 cover all possibilities for relative positions of points X, Y £ F 
such that \ifx — { fy\ < &■ For tpo sufficiently large, it follows that F~ l is a 
Bilipschitz map, for all points X, Y £ F such that \px — <Py\ < 0. 

Case D. Assume points X, Y £ F and \tpx — ¥y\ > &■ Without loss of 
generality we can assume that tpx < <fy- Let us denote with X' = F~ 1 (X) 
and Y' = F -1 (Y). 

First, we consider the case when \(fx~ ( £>y\ < 2ti+9. From (JHJ) we compute 



d{X'X) > d(X f ,0)-d(Y',0) = f{<p x )-f{<pY) 

= f(fx) - f(<Px + (<Py - ipx)) > a'Vx™' 1 - 



(9) 



Next, we consider the case when \ipx — ipy\ > 2n + 9. We extend line OY' 
such that it intersects spiral F' at point Y" £ F' such that ipy" £ [9, 27r + 9) 
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(see Figure E]). It is easy to see that point Y" is well defined. We see that 
f(Y') < f(Y"), because / is radially decreasing, see ([T]). 
Using OH]), we compute 

d(X>,Y>) > d(X',0)-d(Y',0) = f(X')-f(Y') (10) 

> f(X>) - f(Y") = f(tp x ) - f{ipy„) 

From triangle inequality we know that 

d(X',Y') < d(X',X) + d(X,Y) + d(Y,Y'), (11) 

d{X,Y) < d(X,X') + d(X',Y') + d(Y',Y). 

Now, we define (p n := (p(t n ), Vn G No- Notice that from (jSJ) follows cp n ~ t n 
as n — > oo. 

It is easy to see that t 2 k+2 — hk+i < 2n for every k G N . (On the contrary, 
T would not be a spiral, see ([I]), because it would have a self intersection.) 
Using (jHJ) we see that y?2fc+2 — y?2fc+i < 2tt + 1 for every fc 6 No. (For the 
upper bound we could take any number larger than 2ir.) 

Now, assume Z' G r" and Z := F(Z'). There exists kz G No such that 
<fz e b2fc z ,^2fc z +2]- If Vz e [^2fc z ,^2fe z +i] then Z G f |[< 2fcz ,i 2fez+1 ], so 

d(Z, Z') = 0. 

If <Pz e b2fc z +i,^2fe z +2] then Z 6 f |[t 2fcz+1 ,i 2fez+2 ], so (gj(iii) gives 
d(Z,Z') < sup {>(£)- r(£ 2fez+1 )} 

*S[t 2 J; z+1 ,t 2fcz+2 ] 

= ^ fr osc r(t) = o(t 2 ^ + \)=o(^ z - + 1 1 )=o(^- 1 ), 

tt[t2fe z +l,l2fc z +2j 

for y9o sufficiently large. So, for every Z' G V and Z := F(Z') is 

0<rf(Z,Z / ) = o(^ a - 1 ), (12) 

for (/)o sufficiently large. From (TTTj) it follows 
d(X', Y") - d(X, X') - d{Y, Y') < d(X, Y) < d{X', Y') + d(X, X') + d{Y, Y'). 
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Using (|T2|) we get 

d(X', Y') - \o fo^" 1 ) | - \o {ipy^ 1 ) | < d(X, Y), 
d(X, Y) < d(X>, Y') + \o ( Vx a - 1 ) | + \o (<p Y a - 1 ) |, 

for ipo sufficiently large, which is equivalent to 

{ - pT ^ ^d(X',Y') < d(X,Y) <(l+ l °J! Pjrl \ ' ) fHX'.Y') 



d{X',Y') I K ' ' ~ v ' ' ~ \ d(X',Y') 

for (p sufficiently large. (Notice that as ipx < fy, by abusing notation we 
have o ^ipy al ) Q ° ("^x*" 1 )' f° r V^o sufficiently large.) Respecting the lower 
bound on d{X', Y'), p and (|TU|) . we get 

(1 _ \^-)d(X', Y>) < d(X, Y) < (1 + ^)d(X', r), 
a/ a 

for <y9 sufficiently large, which means there exist constants K 3 G (0, 1) and 
K4 > 1 such that 

K 3 d(X',Y') < d(X,Y) < K 4 d(X',Y'), 

for tpo sufficiently large. So finally, for ipo sufficiently large it follows that F~ x 
is a Bilipschitz map on T for points X, Y e T such that |y9x — ^y\ > #• 

From cases A to D it follows that F~ x is a Bilipschitz map on spiral T, 
for <£>o sufficiently large. On the other hand, using (JH)(ii) we see that the 
number of smooth pieces of function f((p) on any bounded interval is finite 
and using (f2k+2 — f2k+i < 2vr + 1, for every fc G N , we see that spiral T, for 
<Po sufficiently large, satisfies all assumptions of Theorem [2j Using Lemma [1] 
and as F is a Bilipschitz map on r", for ip sufficiently large, we finally get 

dim B r' = dim B T 



1 + a 



a 

Finally, using Theorem [3] and Lemma [1] we provide proof of Theorem [1] 
about the box dimension of a spiral generated by generalized Bessel func- 
tions. This proof consists of checking out the conditions of Theorem [3j The 
following lemma makes this checking easier. 

18 



Lemma 2. Let a G (0, 1) and 



r{t)=p(t)\ll + ?^f(t)+g(t), te[t ,oo), t >0, 

where p(t) ~i t~ a as t — » oo, |/(t)| < 1 for all t G [to, oo) and g(t) G 
O (t~ 2 ) as t — > oo. Let C G R and assume that t(<p) = tp + C + 0((/? _1 ) as 
(/9 — t- oo. Let Ay) > 1. 

Then there exists constant k > 0, independent of (p and /Sip, such that 
for all if sufficiently large it holds 

r(t(ip)) - r(t(ip + Aip)) > k<p- a -\l + O^- 1 )). (13) 

The proof of Lemma [2] is omitted, because it is very similar to the proof of 
[3|, Lemma 3]. 

Remark 3. To the contrary to constant k, the constant hidden in term 
O (<^ _1 ) in ffT3|) is not independent of Atp. This is the reason we must pre- 
scribe upper bound on term Aip, in parts of the proof, where we will use 
Lemma EJ 

Throughout the proof of Theorem [1] we will use Hankel's asymptotic ex- 
pansions of Bessel functions J v (t) and Y v (t) for large t, see p, p. 120], as 
follows: 

Mt) = ^Y[Pu(t)co8x-Qv(t)amx], (14) 

Y -(t) = (^) [Pu(t) sin X + Qu(t) cos X ], (15) 

where v G R and x = t — (|z/ + ~) it. P u {t) and Q v {t) are given by 



k=0 

N 



(2ty 



<W*) = Ily^-^^ + oit-^), (17) 
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as t — > oo, which are asymptotic expansions to N terms, where we introduce 
the notation 

(-l) k 1 1 

(4z/ 2 - l)(4z/ 2 - 3 2 ) ■ ■ • (4z/ 2 - (2k - l) 2 ) 



2 2k k\ 
(z/,0) = 1. 

We will also use asymptotic expansions of derivatives 

d ( 2 \ 3 

— J„(*) = I— J [-#„(£) sin* -£„(*) cos*] 



^F,(t) = Q^Y[ik(t) cos* -£,(*) sin*], (19) 

where ^(t) and S„(t) are given by 



N ' ,2 , 1CI.2 



^^-LH)^nF(^F+°( i )> (2°) 



fc=0 

^ ■ ;> i/o/. m2 



as £ — )• oo. 

Notice, that because J^(t) and l^(t) are analytic functions for t > 0, 
see |6|, it is easily verified that -P„(t), Qvif), R v (t) and S^t) are analytic 
functions for £ > 0. This is a sufficient condition for asymptotic power series 
of these functions to be differentiable term by term, see |l|, 1.6]. This property 
will enable us to easily compute expressions containing higher derivatives of 
those functions. 

Proof of Theorem [3 

First, we assume that x(t) := J u ^(t), for every t G [r , oo). The proof is 
presented in a series of steps: 

Step 1. (The box dimension is invariant with respect to mirroring of a spi- 
ral.) We will prove the equivalent claim, that planar curve r" = {(x(t), —x(t)) : 
t G [t , oo)} is a spiral, see (JTJ, defined by r = /(<£>), <p G [<po, oo), near the 
origin, satisfying f((p) — ip~^~, in polar coordinates, near the origin, and 
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dims r" = ^— . It is easy to see that curve T is a mirror image of curve 
r", with respect to the x-axis, hence T is a spiral. Notice that if curve T' 
is a wavy spiral then curve T is also a wavy spiral. Reflecting with respect 
to the x-axis in the plane is an isometric map. As the isometric map is bi- 
Lipschitzian and therefore it preserves box dimensions, see J2|, p. 44], we see 
that dims T = dim^ r" = -^—. 

Step 2. (Checking condition (jSJ).) From 

x(t) = £,„(*) = t^ ./*(*), 

x(t) = J t J ^(t) = ^t—Mt)+t—-Mt) } (22) 



where v = y (^j - ) + v 2 , using ( [141 - 1211 . we compute 

t an f t ) - A(t) - Cgg^j + SrggW) sin * + ($#) ~ fe^(g) cos * 
x(£) P ? (t)cosx - (5p(t)sinx 

(23) 
where we take x = t — \\u + |) 7r. 

Using (12"31) and Lemma [21 since function <£>(£) is continuous by the defini- 
tion, there exists fceZ such that 

^(t) = (x + kn) + 0(t _1 ) as £ ->■ oo. 
For this fc, we define </>o := To + (for — (|^ + t) tm and it holds 

^(t) - O = (t - r ) + 0(t _1 ) as t ->■ oo. (24) 

Define y9 : = v(*o) an d notice that generally O is n °t equal 4> . 

Step 3. (Checking condition (JB]).) From ( 1241 it follows that ip ~ £ as t — > 
oo, and from ( 1221) . using ( TT4H2TJ) we get 



r 2 (£) = (x(t)) 2 + (-x(t)) 2 = -r 2+ ^ + O (r 3+fi ) as £ ^ oo, (25) 

which implies 

f(ip(t)) = r(t) ~ £ 2~ ~ yj 2~ as £ — » oo. (26) 

Notice that from ( |22l) and ( |25|) it follows that function r(£) is of class C 1 . 
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Step 4- (Checking condition (J7|).) By differentiating (|2"3"j) we obtain 

dip.. 9 . . x(t) 2 — x(t)x(t) ,„_ s 

7 | w = orfpW _U_LLLi. (27) 

Using (T2"3"j) again, we have 



cos 2 (o(t) = ^— = ; ^ -. ,„ . (28) 

YK) l + tan 2 y?(£) x(tf + x{tf V ; 

Substituting into ( 127)) and using ffbiH2~T]) we get 

lim -f(t) = lim -^-- — -rrr 1 = L 29 

On the other hand, differentiating (125)) and from (122)) . using ffTH - T2TT) again, 
as — 2 + /j, < 0, we obtain that 

rfr x(t)x(t) + x(t)x(t) f r 3+ ^(-2 + //) sin(x) 2 + g (f 4+M ) 

rft v /x(t) 2 + x(tp _ ~ K0 

~ -r^ ii " 1 sin(x) 2 + (r^^ 2 ) ast^ oo. (30) 

Since f (£) = /'(p) ■ %(t) and since by (J22J we have ^(t) ~ 1 as t -» oo, 
there exists Co, Ci G M, Ci > Co > 0, such that 

\f(<p)\ < C^r^r- 1 < dip-^- 1 as <p -> oo. 

Also, notice that substituting ( 128]) in fl27|) it follows that function (p(t) is 
of class C 1 . 

Step 5. (Checking condition (JSJ).) From (j29|) it follows that there exists 
Ti > r such that -^f(t) > for all t > r l5 so then function </>(£) is strictly 
increasing for all t sufficiently large. As <p(t) is continuous, we conclude that 
for all tp sufficiently large there exists inverse function t = t(<p) of function 
<p = ip(t) and that 

t(<f) = cp + (r - O ) + 0(v9 _1 ) as v? ->■ oo. 

Define 0i := <^(ti) and notice that we can take T\ sufficiently large such that 



?i > <Po- 
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From (1250 . using more terms in the asymptotic expansion, we obtain 
r(t) = ^-p~ v (t)Ml + O (t- 2 )) + (- J|| + O (t-^ sin 2* + O (i~ 2 ) cos 2*, 

where Pv(t) = t~ aP~(t). Using LemmaHwe conclude that for a fixed Aip > 1, 
we have 

f{<p) - f{<p + Ay) = r(t(<p)) - r{t{<p + A<p)) > h^' 1 , (31) 



provided <p is sufficiently large (to use LemmaEl rescaling of r(t) by factor J ~ 

is needed). Moreover, by a careful examination of the proof of Lemma EJ we 
conclude that statement (I3TI) uniformly holds for every Aip from a bounded 
interval, see Remark [31 whose lower bound is greater than 1, provided if is 
sufficiently large. But this means that before using Theorem [31 we will have 
to ensure that 9 from Theorem [3] is larger than 1. It will be sufficient to take 
9 = (tt/3 + l)/2, see Step 9. 

Step 6. (r ; is a spiral near the origin.) Now we can prove that r" is a 
spiral near the origin, that is, /(</?) satisfies condition ([1]) near the origin. 
First, from ([26]) it follows that /(</?) — > as ip — > oo. Second, from (|3"T|) it 
follows that f((p) is radially decreasing for all (p sufficiently large, that is, 
there exists 02 > <Pi such that /|[^ 2j0 o) is radially decreasing. 

S'iep 7. (The box dimension is invariant with respect to taking r and 
(fio sufficiently large.) First, we define r 2 to be such that <^(r 2 ) = 02- Notice 
that r 2 is well-defined and r 2 > T±. From (1231) it follows that r(t) > 0. 
Furthermore, r(t) is strictly positive. Otherwise we would have x(u ) = 
and x(uq) = for some Uq G R and from (J3J) would then follow that x{u) = 0, 
for all u > uq. This is a contradiction with x(t) = J v ,/j,(t) being a nontrivial 
solution of (J3J). As r(t) is strictly positive, there exists constant mi > such 
that for all t€ [tq, t 2 ] it holds 

r(i) > m 1 . (32) 

Notice that 2 > 0i > 0o- From ( 1261) it follows that r(t) — > as t — > oo so 
there exists T3 > r 2 such that for all t G [t^, oo) it holds 

r(t) < mi. (33) 
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We define 3 := f(r 3 ). Notice that we could increase r 3 and 3 to accommo- 
date all requirements in different parts of the proof on t or if being sufficiently 
large. Now, from ( 132]) and (133]) we conclude that 

r / l[n ) ,. 2 ]n r, l^,oo ) = 0. (34) 

As f\u 2! oo) is radially decreasing and function <p(t) is strictly increasing for 
all t G [r 2 , oo), it follows that r'|( T2j00 ) does not have self intersections, so 

r / | [ r 2 ^n r 'l(-3,oo) = 0. (35) 

Finally, from ( 134"]) and ( 135]) we have r'|[ T0iT3 ] flT'l^oo) = 0- Now, we can 
apply Lemma [T]on curve T'. 

Using Lemma [1] we see that without loss of generality we can assume 
that r and <J)q, in the assumptions of the theorem, are sufficiently large. 
Informally, we can always remove any rectifiable part from the beginning of 
T', without changing the box dimension of T'. 

Step 8. (Checking waviness condition (j3|).) By factoring ( 130]) . again from 
(T22D, using flUHHlh we get 

dr [2 

■ (2Pv(t) + 4{t) + cat® (2qL(t) - pg(t))) ^, (36) 

for every \ ¥" k^i k E 'L (division by factor sinx), where pn{t) = £ _ 2P~(£) 
and q„(t) = t~aQv(t). 

Now, further using ([TBI) and (fTTj) . we get 

rtt) + &(*) § 



«#) - J>k(*) 8 - 6/x + ^ 2 + 4z/ 2 

V^)H-^(t) 8(-2 + yU ) 



2^(t) - p£(t) (-4 + /i) (4 - 4/x + /x 2 + 4z/ 2 

Analogously, we see that the derivatives 
d ( Pu(t) + q' v {t) 



t as t — )■ oo, (37) 

£ as t — t- oo.(38) 



* V «/(*) - P'uit) J 8 - 6/i + /i 2 + 4z/ 2 

d / 2pL(t) + 4(t)\ 8(-2 + /i) 



as £ — > oo, (39) 



dt V Hi*) - p£C0 J (- 4 + /i) (4 - 4^ + /i 2 + 4z/ 2 ) 

as £ -> oo, (40) 
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are bounded for t sufficiently large. 

By Lemma [H] and using (ISTHIUI) (possibly with the shift by a constant, 
from variable t to Xi which does not change asymptotic behavior), there 
exists ko G No such that equations 

rtt>+m cot - = _w)+rf(«) (41) 

have unique solutions \2k and X2k-i, respectively, in intervals ((k — 1 + 
k )ii, {k + k )ir), for each A; G N . 

Notice that in (T5TT) and (1551) the right hand side coefficients multiplying 
variable t are strictly less than zero, for every /i G (0,2) and !/6R, So, 
from fl37|) and fl38l) follows that the right hand sides in (J4~T1) tend to minus 
infinity as t — > oo. From the shape of function cot \ it follows that for every 
5 > there exists fco(5) such that solutions \2k and X2fe-i he in intervals 
{{k + fco(^)) 71 " — $, (k + ko(5))ii), for each k G No. Here, we have possibly 
increased k from Lemma [8] and re-indexed sequences \2k and X2k-i- 

Also, notice that if v ^ then the right hand side coefficient in (|37|) is 
strictly less than the right hand side coefficient in ( 138|) . So from (l37j) and (1381) 
follows that if v ^ then x 2 fc 7^ X2fc-i f° r k and analogously £ sufficiently 
large. Contrary, for v = the right hand side coefficients are equal. That 
case is treated in completely different manner and is discussed later at the 
end of the proof. 

In terms of variables t 2 k and t 2 k-i we see that for v 7^ and every 5 > 
there exists k' Q (S) G No (possibly again increased from the value of k (5), to 
accommodate k and t being sufficiently large for X2k 7^ X2k-i and sequences 
re-indexed) such that the equations (J4TI) have unique and different solutions 
t 2 k and t 2 k-i, respectively, in intervals 

(k + k' (5))n + (h + ±) vr - 5, (* + k' (5))n + Qi7 + ±) tt) , (42) 

for each fc G No- 

Without loss of generality we can take £2^-1 < hk- Notice that t 2 k — 
t 2 k-i < 5. For this proof it will be sufficient to take S = 7r/3. It is easy 
to see from ( 136|) that %{t) is positive between these solutions on intervals 
(t 2 k~i,t 2 k), for t sufficiently large (equivalently taking k' (S) to be sufficiently 
large). 

As %(t) > for all t sufficiently large, from § (t) = f'(<p) ■ %(t) it 
follows that f'(<p) > on set UfcLi(v?2fc-i, 02k) where tp 2k -i ■= (pfok-i) and 
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02k '■= f(t2k)- This implies that function f(ip) is increasing for some ip, so 
we can not apply Theorem [2] directly. 

Notice that if t G \JT=o(^k-i,hk) then r'(t) > and if t G Ufc=o(*2fc,*2fe+i) 
then r'(t) < 0. 

We would like to prove that for every k G No there exists unique t 2k G 
(hkihk+i) such that r(t 2k ) = r(t 2k -i) and t 2k - t 2k -\ < n/3 (where we will 
take k' (S) from (JI2"]) to be sufficiently large). As r(t 2k ) > r(t 2 k-i), and as 
function r(t) is continuous and strictly decreasing on interval (t 2k ,t 2k+ i), it 
follows that, if such t 2k exists then it is necessary unique, so we only need to 
prove the existence. 

For every k G N we take t 2k := t 2 k-i + tt/3. Notice that t 2k G (t 2k , t 2k+ i), 
because from (1121) follows that £2*1+1 — hk-i > vr — 5 = 27r/3 and t 2k — 
hk-i < 8 = 7r/3. Define y?2fc := <f{t2k) and take <£>2fe-i as defined before. 
Using (j2"3|) . we can take t or equivalently fc (5) sufficiently large, such that 
(tt/3 + l)/2 < (p 2k — ip 2k -\ < 2 for every k G No- (The exact value of the 
upper bound is not important. We just take some value larger than 7r/3. For 
lower bound, it is only important that it is larger than 1 and lower than 7r/3, 
so we take the mean value between these two.) 

Now, using Lemma[2]and Remark[3j analogously as in Step 5, we compute 

r(t 2 k-i) ~ r(t 2 k) = r(t(<p 2k -i)) - r(t((p 2k )) 

= r(t(ip 2k -i)) - r(t(v9 2fc _i + {(p 2 k - V?2fc-i))) 

> C ^~ 2k ~\~ X > ° 5 

for some C 2 > 0, provided cp or equivalently k' (5) is sufficiently large. From 
this follows r{t 2k ) < r(t 2 fc-i), and as function r(t) is of class C 1 , strictly 
decreasing on interval (t 2k ,i 2k ) and r(t 2k ) > r(t 2k -i), we see that there exist 
hk G {i 2 k,i 2 k) such that r(t 2k ) = r(t 2k -i) and obviously t 2k - t 2k -i < n/3. 

Using t 2k+ i — t 2k -\ > 27r/3, follows that t 2k+ i — t 2k > 27r/3 — n/3 = %/3. 

We established that for every k G No holds t 2k +i > t 2k > t 2 k-x- Notice 
that r'(t ) < and that sequence (t n ), n G No, is the same as the sequence 
from Definition (TJ defined for function r(t). 

As t 2k+ i — t 2k _i > 27r/3 for every k G N , we conclude that t n — )■ 00 as 
n — )• 00, which means that sequence {t n ) satisfies condition (jl])(i). 

As t 2k+ i — t 2k > 7r/3 for every k G N , by taking e = tt/3, we see that 
sequence (t n ) satisfies condition (J4j) (ii) . 

Using ( |30j) we conclude that there exist C3, C4 G 1R, C4 > C3 > 0, such 
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that 

*2fc+2 

osc r(t) = r(i 2 k+2) - r(t 2 k+i) = / r'(t) dt 

*2fc+l 

2 — fl ry ^ 2 — fl 

l 2k+l — °4''2fc+2 ' 

*S[t 2 fc+l,t2fc+2] 



< S ■ sup r'(t) < C 3 t> < C 4 t 2 ~ : 



for every k G N , which means that sequence (t n ) satisfies condition (j4])(iii). 
Finally, we conclude that, for v ^ 0, sequence (t n ) satisfies waviness 
condition (j4]), so r(t) is a wavy function and T' is a wavy spiral near the 
origin. 

Step 9. (Final conclusion.) From the previous steps we see directly that 
for v ^ 0, all of the assumptions of Theorem [3] are fulfilled. We take e' = 
(tt/3 + 1)/2 < e, 9 = min{e', vr} = (tt/3 + 1)/2 and a = ^. Using Theorem[3] 
we prove that 

dim B r' = j^-. (43) 

S'tep 10. (Degenerate case v — 0.) Using known properties and transfor- 
mations of Bessel functions, only in the case of v = 0, expression ( 13 Up for 
function r(t), from ( )22|) . can be simplified in the form 

rj r t^ 4 ( II -21 / \ 2 

-(t) = ^ J (^^(t) + („ - 1 - !„ - 1|) J^,(t)) , 

which is less than or equal to for every t > r . This means that function 
/(<£>) is nonincreasing for (^ sufficiently large, so in this special case, we can 
use Theorem [2] to prove claim ( 1431 . 

The proof for x(t) = Y V4l is analogous and will be omitted. 

□ 



Appendix A. Auxiliary results 

Proposition 1. Let XOY be an isosceles triangle such that R = d(X, O) = 
d(Y,0). Let 6 = Z(XOY) < vr. Then 

6R<-d(X,Y). (A.l) 
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Proof. Using elementary trigonometry we see that d(X, Y) = 2i?sin|. 
Now, it is easy to see that 

n 

6<irsm- } for all G [0,7r], 

which proves (1A.1J) . D 

Lemma 3. (Connection between ip and t) Let t > and v G R. Assume 
if : [to, oo) — > R is a continuous function and let 

, s sin(x) (1 + (r 2 )) + cos( X )0 (r 1 ) 

ttm<p(t) = — 7^- -A ;; , /V ., as£->oo, (A.2) 

yw cos(x)(l + 0(t- 2 ))+sm( X )0(t- 1 ) 

where x = ^ — I — I — ) 7T. T/ien t/iere exists unique k G Z sitc/i t/iot 

yj(t) = fo + A;vr) + O (r 1 ) as £ -> oo. (A.3) 

Proof. Let x is such that | sinxl < ^. Then follows | cosxl > % arid 
| sinxl < | cosxl < 1. From this follows that the denominator of ( 1A.2p equals 
cos(x) (1 + (t -1 )) so expression f|A.2[) becomes 

sin( X )(l + o(t- 2 )) + o(r 1 ) 

cos( X )(l + 0(H)) ast^oo. 

As lower and upper bounds on the denominator are strictly larger than zero 
and strictly lower then infinity, respectively, we write 

g±g + o(r') = o(f) + M,i(i + n) „«-«,. 

Using Lemma H] and Lemma we compute 

tan <p(t) = tan (x + O (t~ l ) ) , (A.4) 

for | sinxl < ^ and t sufficiently large. 

Now let x is such that | sin xl > 2 • Then follows | cosxl < 2 anc ^ 
1 > I sin % I > | cosxl- Analogously as before, we calculate the reciprocal 
expression of expression (|A.2[) to be equal to 

cot <p(t) - cos(x)(l + 0(r 2 )) + sin( X )0(r 1 ) 



sin( X )(l + 0(r 2 )) + cos(x)0(t- 1 ) 

o (r 1 ) + cot(x) (1 + (r 1 )) = cot ( x + o (r 1 )) ,(a.5) 
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for | sinx| > ^ and t sufficiently large. 

Finally, as function (pit) is a continuous function, from (1A.4J) and (IA.5I) 



it is easy to see that there exists unique fcGZ such that it holds (IA.3I) . 

D 
Lemma 4. For every f : [to, oo) — > WL, t > 0, such that lim f(t) = there 

t— >oo 

exist Hi : [to, oo) —¥ K suc/i t/iat y.j(t) = 0(f(t)) as t — > oo, i = 1, 2, and /io/d 

(i) (1 + /(£)) tant = tan(t + yi(i)) /or £ sufficiently large, 
(ii) (1 + /(£)) cott = cot(£ + y 2 (t)) for t sufficiently large. 

Proof. 

(i) For a fixed function f{t) we pointwise show that there exists function 
yi(t) for t sufficiently large. From the definition of the limit it follows 
that there exists t\ > t such that \f(t)\ < |, Vi > t±. 
Now, let t > t\. If /(£) > from Lemma [5] follows that there exists 
y x (t) € R such that \yi(t)\ < f • /(*) and it holds (1 + /(i))tan£ = 
tan(t + yi(t)). If f{t) < 0, as |/(t)| < \, from Lemma [6] follows that 
there exists value y\{t) G M. such that |yi(t)| < n ■ (—fit)) and it holds 
(1 + /(£))tan£ = tan(t + yi(t)). Finally, if f(t) = the statement 
trivially holds for y\(t) = 0. 

So, for constant C = ir, hold \yi(t)\ < C\f(t)\ and (1 + /(t))tant = 
tan(t + yiit)) for t sufficiently large. 

(ii) Define f\it) = fnhk ■ ft * s eas Y to see that lim fAt) = hence, us- 

ing (i), we compute 

(1 + /(£)) COt t = ryrr = ; rTT = COt it + y 2 it) ) , 

v JK " (l + /i(t))tant tan(t + 3/ 2 (*)) 

for £ sufficiently large. 

□ 

Lemma 5. Let a > 0. T/ien /or even/ 1 6 M there exists j/6l snc/j t/iat 
||/| < | • a and it holds (1 + a) tant = tan(t + y). 
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Figure A. 7: Geometry in the proof of Lemma [5l 



Proof. Let t G K 
Case 1. lit 



We distinguish four disjoint cases. 
- kit for some fcGZ. We take y = 0, as it holds 

(1 + a) tan/c7r = = tan(/c7r + y). 

= kit + ^ for some fc G Z. We also take ?/ = 0, because it 



7T 



±oo = tan 



7T 



far + - + y 



Case 2. If t 

holds 

(1 + a) tan (kn, 

Case 3. If t G (&7T, &7T + ~) for some k G Z. We define to '■= t — krc. 

Look at the graph of function tant on interval (k7r,kiT + |), see figure 
IA.7I Let point -E 1 be on the x-axis corresponding to value t, point / be on the 
y-axis corresponding to value tant and point A have coordinates (t, tant). 
Let point J on the y-axis correspond to (1 + a) tant. 

As function tant is a continuous and strictly increasing function on in- 
terval (kn, kn + |) and as statement (1 + a) tant > tant holds, we conclude 
that there exists unique value y such that it holds (1 + a) tant = tan(t + y). 
Let point G on the x-axis correspond to value t + y. As function tant is 
increasing function, notice that y > and y = d(E,G). Also, let point D 
have coordinates (t + y, tan(t + y)). 

Let point C be the intersection of line JD parallel to the x-axis and the 
tangent on the graph of function tant in point A. Because function tant is 
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convex function on interval (kir, kir + ~), point C is to the right of point D 
and have coordinates (t + y + z, tan(t + y)), where z > and z = d(C, D). 
Let point H be the projection of point C to the x-axis, having coordinates 
(t + y + z,0). 

Next, let point B be the intersection of the line JD parallel to the x-axis 
and the line AE parallel to the y-axis. Notice that point B have coordinates 
(t,tan(t + y)). 

Finally, let point F be the intersection of the tangent on the graph of 
function tant in point A and the x-axis. Point F is to the left of point E 
and have the coordinates (t — x, 0), where x > and x = d(E, F). 

Now look at similar triangles AEF and ABC. From similarity it follows 



d(A,E) d{A,B) 



d(E } F) d(B,C)' 



(A.6) 



From ( 1A.6J) we compute 



x y + z V 



tant (1 + a) tant — tant atant 

Finally we get that y < Kx and as function tant is a convex function on 
interval (kir, kn + |), it follows that x < to < f hence y < | • a. As y > it 
holds \y\ < | ■ a. 

Case 4. If t E (kn + f , (k + 1)tt) for some k e Z. Now, -i G (-(fc + 
l)7r, —(A; + 1)tt + |). Using Case 3, we see that 

(1 + a) tant = —(1 + a) tan(— t) = — tan(— t + y) = tan(t + (—y)) 

and that | — y\ = \y\ < ~ ■ a. 

D 

Lemma 6. Let a G (0, |). Then for every t 6 R there exists y £ M. such 
that \y\ < tt ■ a and it holds (1 — a) tant = tan(t — y). 

Proof. Let teR. We distinguish four disjoint cases. 

Case 1 and Case 2. If t — kw or t — /c7r + |, for some k G Z, then we 
take y = 0. The proof is analogous as in Lemma |5j 

Case 3. If t G (krr, kir + |) for some fc G Z. 

Look at the graph of function tan £ on interval (/c7r, /c7r+|), see Figure IAT71 
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Points A to J are defined analogously as in the proof of Lemma EJ The 
important difference here is that point A has coordinates (t — y,(l — a) tant) 
and point D has coordinates (t, tant). 

From similarity equation (1A.6J) . we compute 



x y + z y 



(1 — a) tant tant — (1 — a) tant a taut 

We get that y < -^x < 2Kx and as x < | we get y < n ■ a. As y > it 
holds |y| < 7r • a. 

Case 4. If t G (/c7r + ~, (A; + l)7r) for some fcGZ. Analogous to the proof 
of Lemma [5l 

D 

Lemma 7. For ever?/ / : R — )■ R there exist yi : R — > R snc/i t/iat ?/j(t) = 
0(f(t)) as t — y oo ; i = 1, 2, and hold 

(i) f(t) + tant = tan(t + j/i(t)) /or a// 1 G R ; 
(ii) f[t) + cot t = cot(t + y 2 (t)) for all t G R. 

Moreover, function y\{t) can be chosen such that \yi(t)\ < \f(t)\ for allt G R, 
i = 1,2. 

Proo/. 

(i) For a fixed function /(t) we will pointwise construct function yi(t), by 
defining the value yi(t) for every t G R. 
Let t G R. We distinguish two disjoint cases. 

Case 1. If t = kit + | for some fceZ. We define |/i(t) := 0, as for any 
value f(t) G R, it holds 

/(t) + tant = ±oo = tan(t) = tan(t + yi(t)). 

Case 2. If t G (kir — |, /c7r + |) for some fceZ. Notice that function 
tan(t) is a continuous and strictly increasing function on interval (A;7r — 
~, k'K + |) and its image is set R. Therefore there exists unique value 
ti G (A;tt - f , &7T + f ) such that tan(t x ) = tant + f(t) G R. We define 
yi(t) := ti — t, so we can write tan(t + yi(t)) = tant + /(t). By the 
mean value theorem we get 

/ (t) = tan(t + y\ (t) ) — tan t = tan' £ • y 1 (t) , 

for some value £ G (fc7r — §, &7T + |). As tan'£ > 1, for every £ G 
(A:tt - f,fc7r+ f) it follows that (^(t)] < \f(t)\. 
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(ii) Analogous as the proof of (i). 



□ 



Lemma 8. Let F G C 1 (0, oo) be such that F(z) ~ az as z — > oo for some 
a < 0. Assume that inf F' > — oo. Then there exists a nonnegative integer 
fco such that for each k > ko the equation cotz = F(z) possesses the unique 
solution in J& = (kn, (k + l)7r). 

The proof is given in [3|, Lemma 5] . 
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